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Preface 


This  report  is  taken  directly  from  U.  S.  Naval  Acadetny 
Report  EW-9-75  entitled  "Analog  Simulation  of  a Wave  Activated 
Turbine  Generator  Buoy."  It  was  authored  by  Prof.  Richard  A. 
Hirsch  of  thd  Mechanical  Engineering  Department  MIPR  No. 
Z-70099-5-54326-A.  It  describes  Prof.  Hirsch's  analytical 
development  of  the  subject  simulation. 
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SUMMARY 


This  report  presents  an  analog  simulation  model  of  a Wave  Activated 
Turbine  Generator  Buoy  which  is  currently  under  study  and  evaluation  by  the 
United  States  Coast  Guard.  The  equations  describing  the  heaving  dynamics  are 
developed  as  well  as  the  equations  describing  the  output  of  the  air  operated 
turbine/generator  subsystem.  The  equations  are  scaled  for  analog  computation 
and  the  computer  patching  diagrams  are  presented. 

Although  only  limited  full  scale  data  was  available,  an  analysis  of  this 
data  indicates  that  the  model  developed  will  probably  require  modification. 
Final  conclusions  as  to  the  modifications  required  must  await  further  full 
scale  data.  For  this  reason,  analog  computer  results  are  not  presented  at  this 


! 


w 


TABLE  OF  CONTENTS 

1 . INTRODUCTION 

2.  WAVE  EXCITATION 

3.  WATER  COLUMN  MOTION  EQUATION 

4.  BUOY  MOTION  EQUATION 

5.  AIR  COLUMN  ANALYSIS 

6.  TURBINE  TORQUE 

7.  TURBINE/GENERATOR  EQUATIONS 

8.  ANALOG  SIMULATION 

9.  VALIDATION  OF  THE  SIMULATION  MODEL 

10.  CONCLUSIONS  AND  RECOMMENDATIONS 

REFERENCES 


i 

I 


Page 

1 

4 

11 

14 

21 

30 

36 

39 

65 

71 

72 


? 

i 


1 . INTRODUCTION 


I 

!i 


The  Wave  Activated  Turbine  Generator  (WATG)  has  been  under  study  and 
evaluation  by  the  U.  S.  Coast  Guard  for  some  time.  The  basic  concept  has  been 
described  and  analyzed  by  McCormick  in  References  1,  2,  and  3.  A typical 
configuration  is  shown  in  Figure  1-1. 

Wave  action  causes  heaving  motion  of  both  the  buoy  and  central  water 
column.  In  general,  these  heaving  motions  are  of  different  amplitudes  and 
phase  thus  inducing  relative  motion  of  the  air  above  the  water  column.  The 
air  is  forced  through  nozzles  in  the  turbine  stator  and  causes  the  turbine 
rotor  to  rotate.  The  turbine  rotor  is  directly  coupled  to  the  rotor  of  a 
generator  producing  an  alternating  current  output.  This  output  is  converted 
to  direct  current  which  is  used  to  charge  the  batteries  which  operate  the 
various  navigation  aids. 

The  prior  analyses  have  been  concerned  with  the  effect  of  various  system 
parameters,  associated  with  the  buoy  and  water  column,on  the  power  available 
to  the  turbine  as  a result  of  given  sinusoidal  wave  inputs  of  known  frequency 
and  amplitude.  These  analyses  all  have  the  following  features: 

. The  system  output  is  determined  as  the  power  available  to  the  turbine. 

. The  power  available  is  computed  from  the  unsteady  motion  of  incompressible 
air  in  a straight  air  passage  with  a single  head  loss. 

. The  air  motion  is  determined  by  the  relative  motion  of  the  buoy  and 
internal  free  surface  and  does  not  affect  these  motions. 

. The  buoy  and  water  column  are  each  treated  as  uncoupled,  viscously 
damped,  linear  oscillators  forced  by  sinusoidal  waves  of  different 
amplitudes  but  the  same  frequency  and  phase. 

. The  input  wave  motion  is  unaffected  by  the  motions  of  the  buoy  and 
water  column. 
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The  results  reported  in  References  1-3  are  useful  for  identifying  significant 
parameters  and  in  predicting  general  trends.  For  example,  Figure  5 of 
Reference  2 indicates  that  more  energy  is  dissipated  during  the  relative  rising 
of  the  water  column  than  during  the  relative  falling.  In  addition.  Reference  3 
points  out  that  the  power  available  is  strongly  dependent  upon  the  values  of 
the  damping  coefficients  estimated.  Full  scale  test  data  is  necessary  in  order 
to  identify  the  actual  energy  dissipation  mechanisms  inherent  in  the  system. 

Such  tests  have  been  initiated  by  the  U.  S.  Coast  Guard  but  data  is  not  yet 
available. 

The  study  reported  herein  was  undertaken  in  anticipation  of  the  fact  that 
the  prior  linear  analysis  may  prove  to  be  inadequate  and  also  to  include  the 
turbine,  generator,  and  air  passage  in  the  simulation.  Specifically,  the 
objectives  were: 

. To  model  the  buoy  and  water  column  as  a coupled  system  including  all 
significant  non-linear  effects. 

. To  model  the  actual  air  passage  of  the  hardware  under  test. 

. To  model  the  turbine  and  generator. 

. To  present  the  model  in  a form  suitable  for  solution  on  an  analog 
computer. 

. To  analyze  any  full  scale  data  that  becomes  available  and  modify  the 
model  if  necessary. 

Chapters  2-8  present  the  detailed  development  of  the  simulation  model. 
Chapter  9 presents  the  analysis  of  available  full  scale  data  and  the  conclusions 
and  recommendations  resulting  from  this  study  are  presented  in  Chapter  10. 
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2.  WAVE  EXCITATION 

2.1  Background 

The  excitation  imposed  by  a seaway  on  a floating  body  is  difficult  to 
determine  except  in  the  case  of  simple  bodies.  The  first  simplifying  assumption 
made  is  known  as  the  Froude  - Kryloff  hypothesis  which  states  that  the  presence 
of  the  body  does  not  affect  the  pressure  field  of  the  surrounding  wave  system. 

This  hypothesis  has  been  found  to  give  good  results  for  floating  bodies  and, 
hence,  should  be  applicable  herein.  The  second  simplifying  assumption  is 
known  as  the  St.  Dennis  - Pierson  hypothesis  which  states  that  the  response  of 
the  body  to  a sum  of  wave  components  is  the  same  as  the  sum  of  its  responses  to 
each  wave  component.  This  hypothesis  has  been  found  to  give  good  results  in 
relatively  low  sea  states  and,  hence,  should  be  applicable  to  the  buoys  under 
study  which  are  generally  located  in  semi -protected  waters.  Based  upon  these 
assumptions,  the  response  of  the  buoy  system  to  random  waves  of  small  amplitude  can 
be  found  by  the  superposition  in  random  phase  of  responses  to  steady  state 
sinusoidal  wave  inputs.  Therefore,  the  present  study  will  concentrate  on 
determining  the  system  response  to  a single,  sinusoidal  wave  pattern  derived 
from  classical,  linearized  wave  theory.  Such  waves  are  known  as  Airy  waves 
and  are  analyzed  in  Reference  4,  and  the  results  are  summarized  herein. 

2.2  Airy  Waves 

The  solution  of  the  linearized  wave  equation,  due  to  Airy,  results  in 
waves  of  sinusoidal  profile  and  is  applicable  to  waves  of  small  amplitude 
occurring  in  deep  water.  The  profile  is  given  by 

n(x,y,t)  = cose  (x,y,t) 

where 

= wave  amplitude  at  the  surface  = H/2 
H = wave  height 
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6(Xf  y,  t)  = phase  angle  = k(x  cosa)  + k(y  sina)  - at 

. . Zir 

k = wave  number  " ~ ^ 

X * wavelength  = 

0)  = wave  frequency 
g = gravitational  acceleration 
a = wave  heading 
a = frequency  of  encounter 
t = time 

T = wave  period  = Zir/w 

The  coordinate  system  is  defined  in  Figure  2-1. 

In  this  study,  the  buoy  has  a vertical  axis  of  symmetry  and  all  incoming 
waves  will  be  assumed  to  have  the  same  heading  and  hence,  a = 0.  Also,  since 
the  buoy  is  not  underway,  a = w,  so 

n (x,  y,  t)  = Tim  cos  (kx  - wt) 

Furthermore,  according  to  the  Airy  theory,  the  water  particle  orbital  motion 
is  attenuated  at  any  depth  z and  this  attenuation  is  given  by 


where 


cosh  k(2  + d) 
sinh  kd 


d = depth  to  the  sea  bottom 

For  the  Airy  theory  to  apply, the  depth,  d,  must  be  at  least  half  the  wavelength. 
Now 


cosh  k(z  + d)  * cosh  kz  cosh  kd  + sinh  kz  sinh  kd 
and 

kd  » 2ir  Y 


and  if 
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then 


kd  > IT 


sinh  kd  = cosh  kd 


within  1/2  so 


= nm  (cosh  kz  + sinh  kz) 


sinh  kz  = \ (e''^  - e-"^) 


cosh  kz  > I (e''^  + 


hz  • "M  * 


and,  finally 


where 


n = ^ cos  (kx  - tot) 


D = distance  below  the  still  water  line 


(2-1) 


Equation  (2-1)  is  the  final  form  of  the  wave  profile  used  in  this  study  and 
is  used  to  determine  the  pressure  field  and  thus  the  forces  exerted  by  the 
wave  on  the  floating  body.  Since  the  wave  profile  is  a function  of  the 
coordinate  x,  the  pressure  on  the  body  has  to  be  integrated  over  the  body  in 
order  to  determine  the  total  force.  An  example  of  this  integration  follows. 
2.3  Integration  Of  The  Wave  Pressure 

Employing  the  Froude  - Kryloff  hypothesis,  the  pressure  on  the  floating 
body  due  to  the  wave  is 


where 


p = fluid  mass  density 
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In  the  present  study  only  vertical  motion  of  the  buoy  - water  column  system  Is 
considered  and  If  the  submerged  elements  are  assumed  to  have  straight,  vertical 
sides,  the  total  upward  force  is  found  by  integrating  p^  over  the  appropriate 
area  in  the  x-y  plane  at  depth  D.  For  example,  for  an  annular  body  with  inside 
radius  r and  outside  radius  R,  with  draft  D,  as  shown  in  Figure  2-2,  the  total 
vertical  force  F^  is: 


X = r cose 


dA  = rdrde 


/ / cos  (krcose  - ut)  rdrde 


This  integration  has  been  demonstrated  in  Reference  1 and  proceeds  as  follows; 


cosx  = ^ (e^^  + e"^^ j 


cos  (krcose  - wt)  - j ^ ^1 krcose  ^ ^iut  ^-i krcose ^ 

and  from  Reference  5,  page  355 


where 


J^(z)  = Bessel  function  of  the  first  kind  of  order  n 


^Ikrcose  o i \ 
e de  » 2it  Jjj  (kr) 
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Ztt 

j -ikrcose 


de  = 2ir  Jp  (-kr) 


but,  from  Reference  5,  page  353 


Jq  (-k'^)  = Jq 


so 


r - n H „-kD 
F,  - pg  7 e 


i\ 

Jir  (e'‘‘“^  + e^“^  ) Jq  (kr)  dr 


e"*^^  (cosut)  f (kr)  J (kr)  d (kr) 
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and  from  Reference  6,  page  177 


so,  finally 


[ 


x""'  (x)  dx  = x""'  (x) 


e"'^^  [RJ^  (kR)  - rJ^  (kr)]  coswt 


(2-2) 


Equation  (2-2)  is  used  in  subsequent  chapters  to  determine  the  wave  pressure 
forces  acting  on  the  buoy  and  water  column. 
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3.  WATER  COLUMN  MOTION  EQUATION 


3.1  Introduction 

At  rest  In  still  water,  the  length  of  the  water  column  in  the  center  pipe 
(Figure  1-1)  is  and  assuming  that  the  pipe  has  a constant  radius  R^,  its 
cross  sectional  area  is  A^  = ttR^,  then  the  initial  mass  of  water  is 


If  the  vertical  motion  is  referenced  to  the  SWL  (still  water  line)  and  is 
measured  positive  upward,  let  be  the  absolute  vertical  displacement  of  the 
water  column  and  Zg  be  the  absolute  vertical  displacement  of  the  buoy  and 
pipe,  then  at  any  time,  the  mass  of  water  in  motion  within  the  pipe  is 

\ Cx  * 

and  its  absolute  acceleration  is  z , so  pA^  (L  + - 20)^  = SUM  OF  VERTICAL 

w w w w b w forces 

(POSITIVE  UP) 
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The  vertical  forces  acting  are: 

p Aw  = wave  pressure  force  at  the  bottom  of  the  water  column 
(acting  up) 

-pgA^z^  = weight  of  water  displaced  from  static  equilibrium  (acting  down) 

■P  = air  pressure  force  at  the  free  surface  of  the  water  column 
(acting  down) 

D - = friction  force  at  the  pipe/water  internal  interface  (opposing  the 
'***  relative  velocity) 

Expressions  for  these  forces  are  developed  herein. 

3.2  Wave  Pressure  Force 

The  wave  pressure  force  acting  at  the  bottom  of  the  center  pipe  opening 
is  directly  derivable  from  Equation  2-2  with 

r = 0 


R»Rw 

0,  (0)  * 0 

D - 
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so 


P A 
^nw  w 


= ^-k(L„  - Zj) 


COSojt 


3.3  Friction  Force 


The  frictional  force  due  to  the  relative  motion  of  the  water  in  the  pipe 
is  assumed  to  be  derivable  from  the  Darcy-Weissbach  equation  for  pressure 
loss  in  a pipe  due  to  fluid  viscosity.  From  Reference  7,  page  417,  this 


pressure  loss  is 


. fpIPL  _ 4L 
2D  ■ D ^0 


where 


f = friction  factor 
p = fluid  mass  density 
1 = pipe  length 
D = pipe  length 

lT  = average  fluid  velocity  (over  the  cross  section) 
relative  to  the  pipe 

Tjj  = shear  stress  on  pipe  wall  (acting  on  the  wetted  area  2ttRL) 


For  the  water  column 


R = R. 


L = ^ ^ ^w  ■ 

■ ^w 


*^wB  " ”^0  " 


“wB  ’ T fl-w  * <'b  - ^w>l<^B  - ^w>l 

where  the  velocity-squared  term  is  written  as  shown  in  order  to  give  the 
force  the  proper  sign  opposing  the  relative  velocity.  For  example,  when 
Zw  > Zg  the  water  is  moving  up  relative  to  the  pipe  and  the  force  on  the 
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water  is  down  or  negative.  For  turbulent  flow  In  a moderately  rough  pipe  the 
value  of  f would  be  0.020. 

3.4  Equation  of  Motion 

Combining  the  results  of  the  previous  sections  the  equation  of  motion 
for  the  water  column  Is 

\ - PaA 


* ¥<'■>.  * ^ ■ 
TtpgHR^^Jl  (kR^) 


^w  - ^wH^V  ^w)| 

e-k(Lw  - Zg) 


(3-1) 


The  air  pressure  force,  - Pjy^»  will  be  developed  further  In  Chapter  5. 


4.  BUOY  MOTION  EQUATION 
4.1  Introduction 


The  buoy  Is  assumed  to  execute  pure  heaving  motion  and  thus  the  equation 
of  motion  is 

mg*z*g  = i vertical  hydrodynamic  forces  - D^g  + P30Ag 
where  D^g  is  the  force  due  to  friction  in  the  center  pipe  (developed  in  Chapter 
3)  and  Pg0Ag  is  the  air  pressure  force  to  be  developed  in  Chapter  5.  The 
hydrodynamic  forces  arise  from  pressures  due  to  the  displacement,  velocity, 
and  acceleration  associated  with  both  the  buoy  motion  and  the  wave  motion. 
Following  the  exposition  of  Reference  4,  these  forces  can  be  identified  as: 

WDF  = wave  displacement  force 
BDF  = buoy  displacement  force 
WAF  = wave  acceleration  force 
BAF  = buoy  acceleration  force 
VF  = velocity  forces  due  to  wave  and  buoy 
Expressions  for  these  forces  are  developed  herein,  based  upon  Reference  4. 

4.2  Wave  Displacement  Force 

The  wave  displacement  force  in  the  vertical  (z)  direction  is  found  from 

where  dv  is  a volume  element  and  the  integration  is  over  the  submerged  volume. 
If  it  is  assumed  that  the  buoy  flotation  body  is  an  annulus  with  straight  sides 
and  inside  radius  R^  and  outside  radius  Rg,  and  that  the  pressure  on  the  center 
pipe  thickness  can  be  neglected,  the  above  integral  reduces  to  that  given 
previously  in  Chapter  2.  The  resulting  force  is 

WDF  = e''^^4  ■ ^B^  ‘ 
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where  Lg  is  the  draft  of  the  buoys  flotation  collar. 

4.3  Buoy  Displacement  Force 

The  buoy  displacement  force  is  that  resulting  from  changes  in  displacement 
from  the  static  equilibrium  position. 

BDF  = -pgAgZg 

4.4  Wave  Acceleration  Force 

The  wave  acceleration  force  is  derived  from 


where  is  the  hydrodynamic  inertia  (added  mass)  coefficient  in  heave.  If  it 
is  assumed  that  is  not  a function  of  (x,  y,  z)  or  that  some  mean  value  is 
used,  then 

where  F^  is  given  in  Chapter  2,  The  final  form  is 

WAF  = y-  e~'^('-B  ■ ^b)  [^ReJilkRe)  - R^J,(kR^)j  cosu)t 

4.5  Buoy  Acceleration  Force 

The  force  due  to  buoy  acceleration  is 

BAF  = - pVC^  Zg 
but  since  ^ ^*'B  " ^B^ 

BAF  * -pA0Cj^  (Lg  - Zg)  Zg 

4.6  Velocity  Forces  due  to  Wave  and  Buoy 

All  of  the  velocity  forces  are  due  to  relative  velocity  between  the  buoy 
and  the  wave.  There  are  three  distinct  types  of  velocity  forces  to  consider; 
friction  force,  pressure  force,  and  radial  dissipation  or  wavemaking  force. 
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4,6.1  Friction  force, 

Fluid  viscosity  gives  rise  to  a shear  force  on  the  external  vertical 
surface  of  the  buoy  which  is  proportional  to  the  wetted  area  and  the  square 
of  the  relative  velocity.  This  force  is  generally  much  smaller  than  the 
other  velocity  forces  and  hence  it  will  be  neglected. 


= 0 


4.6.2  Pressure  force, 


This  force  is  due  to  the  dynamic  pressure  acting  on  the  cross-sectional 
area  normal  to  the  flow.  The  general  form  is 

Op  . I p Cp  A 

The  area  normal  to  the  vertical  flow  is  Ag  and  the  relative  normal  velocity  at 
the  location  of  this  area  is 


From  Chapter  2 


= n - Zg  at  depth  - Lg 


wH  -k  ( Ln  - 2„)  . 


B " B'  sin  (kr  cose  - wt) 


n - ^ 

In  order  to  preserve  the  proper  sign  for  this  force,  it  must  change  sign  when 

the  relative  velocity  changes  sign.  It  is  also  to  be  noted  that  since  the 

buoy  is  only  partially  submerged,  the  proper  drag  coefficient  is  the  forebody 

drag  coefficient  when  the  relative  velocity  is  postive  and  the  base  drag 

coefficient  when  the  relative  velocity  is  negative.  In  view  of  the  foregoing 
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Dp  = 2 P ^B  "9"  '' 


Now,  since  n varies  over  the  area  (it  is  a function  of  x,  y)  an  integration 
over  the  area  must  be  carried  out,  so  that 


7 p Cj.  (sgn  v) 


Rg  2it 


/ /o  (^  - "b)'  = I 
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I = n^rdrde  - Zz 


hfj  ^ jj 


rdrde 


^ 4 ^3 


There  integrals  are  now  evaluated. 
4.6.2. 1 evaluation  of  I3 


R Ztt 
' B c 


I3  -/  / rdrde  “ Aj  - , - r 

0 


4. 6. 2. 2 evaluation  of  I, 


. |H  ^-KIb  - Zj)  j-B 


Rn  2Tr 


sin  (krcose  - wt)  rdrde 


«w  ° 


in  (krcose  - a,t)  » ^ (krcose  - ait)  _ ^-i  (krcose  - o,t^ 


gi krcose  _ e^“^  g-ikrcose 


/\dr  ff 

\ »L 


e gikrcose  _ e^“^  .-ikrcose" 


fe-MLe-V  / 

- -p-  e ^4  4^  (sin  ut)  J*  kr  (kr)  d (kr) 

«w 
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and  finally 


® [RgJ^  (kRg)  - R^Jj  (kR^)]  sinut 

4. 6. 2. 3 evaluation  of  I, 


^ - 


-1^  Rp  2it 


sin^  (krcose  - tot)  rdrde 


but 


stn2  (krcose  - «.t)  » - | >'(krcose  - »t)  ^ ^-2i(krcose  - «t)  . jJ 


1 e-2'“‘  7, krcose  e^'“‘ 

?■ T ® - J 


SO 


I, 


fwH  g-k(Lg  - Zg)  1 


_ 1 


IT  COS  2 (ut 


7 ■ 7 ^'B  ■ -JWP 


J (2kr)  Jq  (2kr)  d (2kr) 


w 


" 1 ^B  ■ fk  ^Vi  ■ Vi  ^2kR„)]  COS  2 „t 


and  finally 


I,  = 


1 " I 


■ ^B^ 


A - - 
''b  k 


RgJj  (2kRg)  - R^Jj  (2kR^)  cos  2 wt 


J cos  2 u)t  j 


The  final  form  of  the  pressure  force  is 

Dp  = ]p  (Ij  - 2Zgl2  + Zg^Ig)  (Cp  sgn  v) 

where 


sgn  V = 


_ ^ e"*'^*'B  " ^B^ 
B 


Vi  "'“b'  - Vi 


sinait  - z 


B 


and  I3  are  as  above. 

4.6.3  Radial  dissipation  force.  Dp 

This  force  is  generated  by  the  energy  carried  away  by  the  waves  created 
by  the  oscillating  buoy.  For  a hull  with  a vertical  axis  of  symmetry  and 
straight  sides,  the  radial  dissipation  generated  by  a section  of  unit  length  is 
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J 


^ (sin2ky)  ‘ (f,  - Zg) 

where  y Is  the  half-width  of  the  section  and 

n * # ■ ^B^  sin  (kx  - a>t) 


The  total  force  is  then 


Dr  " ^ e"^''^‘‘B  ■ ^B^ 


J (sin2ky)  1^-  e‘'^^'-B  " ^B^  sin  (kx  - (ut)  - ZgJ 


where  the  integration  is  over  the  length  in  the  x-di recti on  and,  in  general,  y 
is  a function  of  x.  This  integration  could  not  be  carried  out  exactly  for  the 
case  at  hand  where  the  water  plane  area  is  circular.  The  integral  can  be 
approximated  by  replacing  the  buoy  waterplane  with  an  equivalent  area  rectangle 
of  width  w and  length  L such  that 

Ag 

y = I = Rg  and  L = 

Proceeding  with  this  approximation 

L/2 

Dr  = ^ e’^'^^'-B  ■ ^B^  sin2kRg  f e^'^^^B  “ ^B^  sin(kx  - oat)  - Zg  dx 

-L/2  ^ 

= ^ sin2kRg  I - e‘''^‘-B  " ^B^  j^cos  - u)t)  - cos(|^  + oat^ 


[*  LA  A 

f ,-'‘('■6  - sin  ^ s1n.t  * ijJ 

which  is  the  final  form  of  the  radial  dissipation  force. 

4.7  Equation  Of  Motion 

Combining  the  results  of  the  previous  sections  the  equation  of  motion  for 


the  buoy  is 


p^b'-b  ^ ^ S ( ^ • q ) J ^b  “ • Dwg  ■ 


■(^sih^icrJ  .-'‘"■b  ■ 'e>  23 
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* r * "bV*  <‘^D  ^9"  ''!  * F I,  (C„  sgn  v)  t p^jAj 


* e‘'‘<'-B  ■ ^B> 


1 - C|^k(L0  - Zg)  I coswt 


If  A “1 

e-S^fl-B  - ^B>  sm^kRj  s1n  ^ 


si  nut 


(4-1) 


where 


^ " I 


A„  - 


7T 

B ■ F 


RgJj(kRg)  - R^J^(kR^)jcos2utj  (4-2) 


I2  = - - Zq)  RgJ^(kRB)  - R^,Oi(kR„) 


s i nut 


(4-3) 


\ = ■ ^B^  (^B  ■ I - ^w) 


(4-4) 


sgn  v = 


ttmH  _-k(Lp  -Zp) 
kAg  ^ 


RbJ.IkRb)  - Vi'i-'w’ 


sinut  - z 


(4-5) 


The  air  pressure  force  on  the  buoy,  P^gAg,  is  developed  in  the  next  chapter. 
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5.  AIR  COLUMN  ANALYSIS 


I 


5.1  Introduction 

The  motion  of  the  air  above  the  water  column  free  surface  is  analyzed  in 
order  to  determine  the  mass  flow  delivered  to  the  turbine  and  the  pressure 
forces  acting  on  the  water  column  and  buoy.  The  question  arises  as  to  just 
how  to  treat  this  air  flow;  compressible  or  incompressible,  viscous  or  non- 
viscous,  steady  or  unsteady?  The  assumptions  made  are  that  the  flow  is 
incompressible,  non-viscous,  and  steady.  These  assumptions  are  justified 
by  the  arguments  to  follow. 

The  air  at  the  water/air  interface  has  velocity 


and  this  velocity  will  rarely,  if  ever,  exceed  a few  feet  per  second  for  the 
buoy  systems  and  deployments  of  interest.  Throughout  the  air  passage,  which 
is  shown  in  Figure  5-1,  compressibility  effects  can  certainly  be  neglected. 
The  turbine  nozzle  areas  are  such  that  the  air  velocity  will  increase  by  a 
factor  of  about  150  and  so  even  at  the  nozzle  exits  the  velocity  will  rarely 
exceed  several  hundred  feet  per  second  and  compressibility  effects  can  again 
be  neglected. 

The  fluid  can  be  treated  as  non-viscous  but  in  computing  pressures 
various  head  losses  can  be  included.  With  these  assumptions,  the  Bernoulli 


equation  along  any  streamline  is 


r 3V  V ' P2  Pi 


and  for  the  moment  consider  a straight  passage  from  the  internal  free  surface 
to  the  turbine  nozzles.  Since  the  continuity  equation  applies,  3V/3t  is  not 
a function  of  s,  the  distance  along  the  passage,  so 


21 


V * p.  V 2 p 

V(z^  - z^)  + 4-  + + g(z^  . z^)  = -2_  + _ 


Now,  assuming  the  nozzles  exhaust  to  atmosphere  and  applying  the  continuity 


equation 


= (Vt9)(z,  (sj) 


where 


Now  assume  a steady-state  response  to  a sinusoidal  wave  Input  so  that 

V = An  slnOt 


V = An*  cosnt 


Pi  A2 

— = (An*  cosnt  + g)  L.  + ^ 

P Q d 


where 


An  = magnitude  of  air  velocity 
= length  of  air  passage 

oi 


approximate  values  are 


L,  = 3 feet 

d 


/A 

I " ] 

VO  ■ 


1 = 21400 


“ 3(An*  cosnt  + g)  + 10720{An)*  sln*nt 


For  a typical  air  motion  values  might  be 


i 


Afi  = 4 FPS 
n = IT  RAD/SEC 


so 


Pi 

— 100  + 40  cosSlt  + 160000  sin^nt 

P 


so  that  the  only  time  at  which  the  unsteady  and  gravity  terms  have  an  ap- 
preciable effect  Is  In  the  vicinity  of  Jit  = 0 and  180°,  and  at  these  times 
the  pressure  Is  small  (on  the  order  of  1/4  PSF).  As  far  as  the  pressure 
forces  are  concerned,  then,  these  terms  can  be  neglected.  These  terms  will 
have  an  effect  on  the  time  at  which  the  Intake  valve  opens  or  closes  (Pj  = 0) 
but  the  precise  timing  of  this  event  Is  not  considered  significant.  In  fact, 
the  valve  Is  a simple  flapper  of  very  low  Inertia  and  It  will  be  assumed 
that  when 

> 0 valve  Is  closed 

z^  < 0 valve  Is  opened 

With  these  assumptions,  then,  the  air  velocity  delivered  to  the  turbine  Is 


V = J 

J IT  A 

for  > 0 

(5-1) 

for  z^  < 0 

(5-2) 

The  air  pressure  forces  can  now  be  determined  by  the  application  of  the 
steady  form  of  the  Bernoulli  equation. 

5.2  Intake  Valve  Closed  (z^  > 0) 

In  the  analysis  that  follows  the  symbols  used  are 

p = air  pressure  p = air  density 

V = air  velocity  h = loss  coefficient 
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and  the  subscripts  refer  to  stations  In  the  air  passage  which  are  Identified 
In  Figure  5-1. 


5.2.1.  Station  1 to  2 


P,  V * p^  V " 


r /Ai\' 

Pi  - P.  = -T-  0 + h,,)  {^) 


h = loss  coefficient  for  a sudden  contraction 


5.2.2.  Station  2 to  3 


D V * D V * 


A,V.  = = A,V, 


Pi  - P,  ' -5-  ('  + l>i,)  (j^) 


h = loss  coefficient  for  a convergent  section 


5.2.3.  Station  3 to  1 


There  Is  no  area  change  In  this  section  and  If  friction  at  the  wall 


Is  neglected 


Pa  * Pl 


5.2.4.  Station  1 to  j 


Pi  ''l*  , Pa  ^a* 

-^  + -^(1  +h,j)  --^  + -2- 
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AjVj  . A,V,  = A,V.  . A,V, 

pv,=  r /’A,\> 

pj  = (1 1 h,j) 


h . = nozzle  loss  coefficient 

* J 


5.2.5.  Air  pressures 

Assuming  that  the  nozzles  exhaust  to  atmosphere,  and  since  Vj 

and  A = A the  final  pressures  are 
I w 

.ft)' 1 


rc''' 


Pa  = Pi 


pV 

T-  VST; 
L J 


.tj  ^ VsT 


5.2.6.  Force  on  water  column 

The  air  pressure  force  on  the  water  column  acts  down  (-)  during  this 


phase  (z^  > 0)  and  Is 


= - Paw  \ - P/w 


^aw^  ~ \ 


(feji  »■" 


This  is  the  air  pressure  force  required  In  Equation  (3-1). 


5.2.7.  Force  on  the  buoy 

The  air  pressure  force  on  the  buoy  acts  up  (+)  during  this  phase 


(z^  > 0)  and  Is 
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'"b  = PaB^B  ' P2^^  • ^ Ps^'^s  - + P23^^  - ^s)  + '"m 

where  Pj,  is  the  area-averaged  pressure  in  the  short  convergent  section  between 
stations  2 and  3 and  is  the  force  due  to  the  turning  of  the  air  stream 

180®  as  it  is  exhausted  from  the  annular  outlet  of  area  A^.  It  can  be 
shown  that 


This  is  the  air  pressure  force  required  in  Equation  (4-1). 


5.3  Intake  Valve  Open  (z^  < 0) 

A similar  analysis  is  made  for  the  intake  flow  when  the  valve  is  open. 
5.3.1.  Station  a to  i 

If  the  velocity  on  the  atmospheric  side  of  the  intake  valve  is  assumed 
to  be  zero 

A.V.  = A V 
i i 11 


PV*  M 

Pi  - Pa r ^^iHy 
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= loss  coefficient  for  the  inlet  port 


5.3.2.  Station  i to  3 


P3  - Pi  = 


5.3.3.  station  3 to  2 


pv, 


P2  - P3  = 


2 -• 


h32  = loss  coefficient  for  a diverging  section 


5.3.4.  Station  2 to  1 


Pi  - P2  = 


PV,= 


'A  Y 

<r)  - ^l) 

2 


= loss  coefficient  for  a sudden  enlargement 


5.3.5.  Air  pressures 


Assuming  that  p,  is  atmospheric,  the  final  pressures  are 

a 


Pi  = 


P2  = 


pz. 


1 + 


^ {^1  ^ K.  m 


pz. 


p,  = - 


* "al 


Pi 


5.3.6.  Force  on  water  column 

The  air  pressure  force  on  the  water  column  acts  up  (+)  during  this 
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1 


phase  (z^  < 0)  and  Is 

fw  = PaA  ' P.*.  ' - K 


(5-5) 


This  is  the  air  pressure  force  required  in  Equation  (3-1). 


5.3.7.  Force  on  the  buoy 

The  air  pressure  force  on  the  buoy  acts  down  (-)  during  this  phase 
(z^  < 0)  and  is 

•"b  = PaB^B  = P2(^  - + Pas^^  " ^3^  Pi^3 


where  Pjs  is  the  area-averaged  pressure  in  the  short  divergent  section  between 
stations  3 and  2.  It  can  be  shown  that 


P23  ~ P3 


so 


PaB^  " ■ ^B 


■2  ^B 


(5-6) 


This  is  the  air  pressure  force  required  in  Equation  (4-1). 
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6.  TURBINE  TORQUE 
6.1  Basic  Theory 

Basic  reaction  turbine  theory  from  Reference  8,  Chapter  5 is  presented 
in  order  to  analyze  the  air  flow  over  the  turbine  blades  so  that  the  torque 
produced  can  be  related  to  turbine  geometry  and  buoy  system  variables.  The 
configuration  of  the  velocity  vectors  is  shown  in  Figure  6-1.  The  symbols 
are  defined  as  follows: 


Vj  = velocity  of  the  air  exiting  the  nozzles 
U = relative  air  velocity  due  to  turbine  rotation 


V-j  = net  entering  air  velocity 
^2  - net  exiting  air  velocity 
= mean  air  velocity 
Oj.  = nozzle  angle 
d-i  = entering  angle  of  flow 
&2  - exiting  angle  of  flow 
= mean  angle  of  attack 

= tangential  force  on  a blade  due  to  lift  and  drag 


From  Reference  8 


h = -7-  “m  - 

Ajj  = blade  area  = span  x chord  = be 


Reference  8 states  that  for  well -designed  blades  the  drag  terms  are  generally 


negligible  so  this  will  be  assumed.  The  lift  coefficient  can  be  approximated 


as 


= 2'fK  sin  Ojjj 

where  K is  the  cascade  coefficient  and  depends  upon  the  outlet  angle,  62, 
and  the  blade  pitch-to-chord  ratio  (s/c).  Substituting 


30 


pV, 


m 


A.  (2TrK)  sin  a„  cos  a 
D m IT 


Now,  from  the  velocity  diagrams  of  Figure  6-1 


sin  = -V2  cos  &2  Vi  cos  ai 
m m ^ ‘ ‘ 


V„  cos  a„  = V2  sin  62  + Vi  sin  ai 

m m ^ ^ ^ 


and,  to  account  for  losses,  let 


V2  = Vi 


h|j  = blade  loss  coefficient 


so  that 


sin  ctjjj  cos  = Vi^  (-h^  cos  62  + cos  aj)  (hj^  sin  &2  ^ sin  aj) 


Also,  from  the  velocity  diagrams 

2 7 2 

« I#  ^ 11^ 


V,  = V/  + U - 2V,U  cos  a. 


J 


V.  cos  a-  - U 

M J W 


COS  aj  = ^ , Sin  ai 


V-  sin  a. 
-J i 


and 


(6-1) 

(6-2) 


= irpKAjjVi  (-hjj  COS  62  + cos  aj)  (hj^  sin  $2  + sin  ai)  (6-3) 

Equations  (6-1),  (6-2),  (6-3)  allow  the  turbine  torque  to  be  found  as  a 
function  of  air  velocity  and  turbine  angular  velocity  since 

T = F,R^ 


U = RjUi 


V = — 2 
J Aj  a 
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where  Rj  is  the  turbine  radius.  From  an  analog  simulation  standpoint,  these 
equations  could  be  solved  but  would  require  seven  multipliers.  Since  these 
equations  are  not  precise,  it  seems  desirable  to  take  an  approach  that  would 
require  fewer  multipliers  but  still  preserve  the  essentials  of  the  torque 
variation. 

6.2  Simplified  Torque  Analysis 

For  an  ideal  impulse  turbine  the  torque  - speed  variation  is  linear  with 
maximum  torque  at  zero  speed  and  zero  torque  at  an  RPM  corresponding  to  blade 
stall.  Reference  8 states  that  a well-designed  reaction  turbine  approximates 
this  behavior.  This  approximation  will  be  employed. 

If  the  turbine  has  n blades,  then  at  zero  RPM  U = 0 and  Equations  (6-1), 

(6-2),  (6-3)  can  be  reduced  to  give  the  total  torque  T^  as 

T = C V ^ 

'o  Vj 


Cj  = nuptCAjjRy  (-hjj  cos  &2  cos  aj)  (hj^  sin  B2  + sin  oj) 


and,  employing  the  linear  approximation 


T = Tq  - mu) 


where  m is  to  be  found. 

From  Equation  (6-3)  the  torque  is  zero  when 


cos  ai  = h.  cos  62  * 


V4  cos  Oi  - U. 


so,  using  Equation  (6-1) 


(Vj^  + - 2VjUj,  cos  Oj)  {h^  cos  cos^  a.  - 2VjUj,  cos  oj  + 


2 2 cos2  a.  - (h^  cos  B2) 

- (2V.  cos  a.)  + V/  — r-  ^ n \z = 0 

0 ' J J 0 j ' “ '"b  ^2' 
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k = 1 - (tan  a.  )Ji  (H -cos  e, 


(6-6) 


Cy  = nirpKAjjR^  (-hjj  cos  pj  + cos  oj)  (hj^  sin  ■*■  sin  o^)  (6-7) 


^ vJ 


for  > 0 


(6-8) 


Vj  = 0 


for  z^  < 0 


(6-9) 


This  set  of  equations  gives  the  input  torque  to  the  turbine/generator  system 
as  a function  of  the  buoy/water  column  motion. 


7.  TURBINE/GENERATOR  EQUATIONS 

7.1  Generator  Characteristics 

The  WATG  employs  a 3-phase  alternator  to  generate  an  a.c.  signal  which 
is  rectified  and  supplies  current  to  charge  a 12-volt  battery.  If  the  battery 
is  fully  charged  a control  circuit  diverts  the  current  to  a resistor  where  it 
is  dissipated.  The  no-load  voltage  and  charging  current  into  a 12-volt  load 
are  shown  in  Figure  7-1  as  a function  of  generator  RPM.  The  generator  is 
directly  coupled  to  the  turbine  so  generator  and  turbine  RPM  are  the  same. 

The  detailed  analysis  of  the  electromechanical  equations  is  quite  complex  but 
does  not  appear  to  be  necessary  for  this  study  since  only  two  questions  need 
be  answered:  First,  given  the  turbine  angular  velocity  what  is  the  current 
generated?  Second,  given  the  input  torque  due  to  air  flow  what  is  the  tur- 
bine angular  velocity?  The  first  question  is  answered  by  using  the  generator 
characteristic  shown  in  Figure  7-1.  For  the  analog  simulation,  this  curve 
is  duplicated  using  a variable  diode  function  generator.  The  second  question 
is  answered  by  analyzing  the  turbine/generator  dynamics  using  a simplified 
equation  of  motion  which  includes  an  electromagnetic  torque  proportional  to 
the  current  generated.  The  equation  of  motion  is  developed  in  the  next  section. 

7.2  Turbine/Generator  Equation  of  Motion 

The  equation  of  motion  is  developed  from  the  relation 

Ju  = sum  of  torques 

J = mass  moment  of  inertia  of  turbine  rotor  and  generator  rotor 
The  torques  applied  consist  of 

T = Torque  due  to  air  flow  over  the  turbine  blades  (given  by 
Equation  (6-4)) 

Tq  = torque  due  to  friction  and  windage  losses 

T^  » torque  due  to  current  generated 
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The  torque  due  to  friction  and  windage  losses  is  assumed  to  be  directly  pro- 
portional to  angular  velocity  and  opposed  to  it.  The  constant  of  proportion- 
ality is  designated  B and  is  equivalent  to  a viscous  damping  coefficient. 

The  simplest  representation  of  the  electrical  torque  is  that  it  is  directly 
proportional  to  the  current  generated  and  the  constant  of  proportionality  is 
kg  the  slope  of  the  no-load  voltage  curve.  With  this  assumption,  the  equa- 
tion of  motion  is 

Jj)  = T - Bo)  - kgi  (7-1 ) 

where 

T is  given  by  Equation  (6-4) 

B = viscous  torque  coefficient  (FT-LBS/RAD/SEC) 

kg  = generator  no-load  voltage  constant  (FT-LBS/AMP) 

i = f(a))  given  by  Figure  7-1  (7-2) 
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8.  ANALOG  SIMULATION 


i 

I 

I 


i 


8.1  Introduction 

The  equations  necessary  to  simulate  the  WATG  have  been 
preceeding  chapters  and  are  identified  as  follows: 


ELEMENT 

Water  Column  Dynamics 
Buoy  Dynamics 
Air  Pressure  Forces 
Turbine  Torque 
Turbine/Generator  Dynamics 


developed  in  the 

EQUATION  NUMBERS 
(3-1) 

(4-1)  through  (4-5) 

(5-1)  through  (5-6) 

(6-4)  through  (6-7) 

(7-1)  and  (7-2) 


The  scaled  computer  equations  and  the  analog  patching  diagrams  are  developed 
in  this  chapter.  The  symbols  used  in  the  diagrams  are  shown  in  Table  8-1. 

8.2  Wave  Forcing  Functions 

An  examination  of  equations  (3-1)  and  (4-1)  through  (4-5)  shows  that  wave 
forcing  functions  proportional  to  sinOt,  cosJlt,  and  cos2fit  are  required.  The 
generation  of  unit  amplitude  sinusoidal  functions  is  a standard  analog  tech- 
nique and  the  required  patching  diagrams  are  shown  in  Figure  8-1.  The  charac- 
teristics of  the  input  wave  are  shown  in  Table  8-2  where 

T = wave  period 

n = 2Tr/T  = wave  frequency 
oT* 

A = = wavelength 

k = 2Tr/X  = wave  number 

Cti  = [kRBJi(kRB)  - kRwJ,(kRw)] 

and  j is  the  amplitude  of  the  wave  at  a depth  equal  to  the  buoy  draft  Lg. 
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Table  8-1.  ANALOG  SYMBOLS 


Table  8-1.  CONTINUED 


ELEMENT 


Relay 

Comparator 


Zero 

Limiter 


Interconnection 


Output 


SYMBOL 


Zl 

Z2- 
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-<i) 


■0 


EQUATION 


If 


If 


X > 0 

y = z, 

X < 0 

y = z. 


If 


X > 0 
y = -X 


If  X < 0 
y « 0 


ShoMS  connections  among 
various  circuits 


ShoMS  recorder  channel 
number 


Table  8-2.  WAVE  EXCITATION  PARAMETERS 


Period  T 
(sec) 

Frequency  £1 
(rad/sec) 

Wavelength  X 
(ft) 

Wave  Number  k 
(ft-*) 

Amplitude  Coefficient 

Cn 

1.0 

6.283 

5.125 

1.226 

-0.0014 

1.2 

5.236 

7.380 

0.8514 

0.0016 

1.4 

4.488 

10.04 

0.6255 

0.0295 

1.6 

3.927 

13.12 

0.4789 

0.0857 

1.8 

3.491 

16.60 

0.3784 

0.1592 

2.0 

3.142 

20.50 

0.3065 

0.2384 

2.2 

2.856 

24.80 

0.2533 

0.3156 

2.4 

2.618 

29.52 

0.2128 

0.3869 

2.6 

2.417 

34.64 

0.1814 

0.4510 

2.8 

2.244 

40.18 

0.1564 

0.5076 

3.0 

2.094 

46.12 

0.1362 

0.5573 

3.2 

1.963 

52.48 

0.1197 

0.6008 

3.4 

1.848 

59.24 

0.1060 

0.6387 

3.6 

1.745 

66.42 

0.0946 

0.6720 

3.8 

1.653 

74.00 

0.0849 

0.7012 

4.0 

1.571 

82.00 

0.0766 

0.7269 

4.2 

1.496 

90.40 

0.0695 

0.7495 

4.4 

1.428 

99.22 

0.0633 

0.7696 

4.6 

1.366 

108.4 

0.0579 

0.7875 

4.8 

1.309 

118.1 

0.0532 

0.8034 

5.0 

1.257 

128.1 

0.0490 

0.8177 

5.2 

1.208 

138.6 

0.0453 

0.8305 

5.4 

1.164 

149.4 

0.0420 

0.8420 

5.6 

1.122 

160.7 

0.0391 

0.8524 

5.8 

1.083 

172.4 

0.0364 

0.8619 

6.0 

1.047 

184.5 

0.0340 

0.8705 
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The  values  used  are  for  an  8 x 26  buoy 

Rg  = 4 feet 

R,  = 1 foot 
w 

Ag  = TT(Rg2  . R^2)  = 15^  ^^2 
Lg  = 4 feet 


8.3  Water  Column  Dynamics 

For  simulation  purposes  equation  (3-1)  is  simplified  by  assuming  that  the 

vertical  motions  of  the  buoy  and  water  column  are  small  compared  to  the  water 

column  length.  Dividing  through  by  pA  L and  letting  the  scaling  factors  be 

w w 

denoted  by 


A = acceleration  scale  factor 
V = velocity  scale  factor 
6 = displacement  scale  factor 
The  scaled  equation  is 


^w 

X 


•^aw 

pL/ 


+ H , cosnt 
w 


where 


H 

■ I 


e‘*^*-w 


_ fv^ 

-4R7 


R^  = 1 foot 
w 

f = 0.020 
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AH  -kL^, 

-j^=  32.2  J^(k)  (sec-*) 


= 0.005 


Since  the  water  column  length  Is  to  be  a parameter  In  the  simulation.  Table  8-3 
presents  these  coefficients  for  various  values  of  L^.  The  air  pressure  force 
term  is  given  in  Section  8.5.  The  analog  patching  diagram  is  shown  in 
Figure  8-2. 

8. 4 Buoy  Dynami cs 

For  simulation  purposes,  equations  (4-1)  through  (4-5)  are  simplified  by 
assuming  the  vertical  motions  of  the  buoy  and  water  column  are  small  compared 
to  the  buoy  length  Lg.  Dividing  through  by  pAgLg(l  + Cj^)  and  using  the  pre- 
viously defined  scaling  factors,  the  buoy  equation  of  motion  is 


Cnsgn  ^ 


7 

7 

/ 

B - r 

■ ■ ''BVj 

^B 

T_ 

+ 

Cbv^  (2Hgi  sinat 

■"b 

^w 

in  1 

.T‘ 

V 

Jl'  V ■ V ■ 

+ Cgj^  ffg  cosnt  - Cgj^  Rg  sin«t 


i^”B2  ( 


Hb2  - Hg^  cos2nt 


)(coS9n|) 


pLbA(1+C 


where 


Y • Rg^  sinflt  + Y 
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Table  8-3.  COEFFICIENTS  FOR  WATER  COLUMN  EQUATION  OF  MOTION 


r 


k = 30  ft 
w 

•c 

o" 

c 

<o 

1.073 

1.073 

a’ 

h 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0003 

0.0009 

0.0023 

0.0049 

0.0090 

0.0148 

0.0222 

0.0314 

0.0420 

0.0538 

0.0667 

0.0802 

0.0943 

0.1087 

0.1232 

0.1377 

0.1520 

0.1660 

0.1798 

0.1932 

LT) 

C\J 

II 

S 

*0 

2 

o 

<c 

«o 

1.288 

1.288 

'5l= 

OOOOOrOr-r-0%0>COCNJ<'^OOOCSJOJC\JCNJCy>CMOrOO>00 

oooooor— covocMi— cMtoor^^rocvji— ocor^mcvjoo*:*' 
ooooooooor— cs4ro<i-vor^o\r-coi^>r*^ooocM^ir>r^ 

OOOOOOOOOOOOOOOOi—  1—  r—  CsJCsJCSJCVJCSJ 

odoooooodooooooooodooooooo 

L = 20  ft 
w 

tc 

2 

o 

c 

►1 

K) 

019*1 

019*1 

'5* 

H 

OOOO^r^ororocvjr^cn^cvjcNjr^'si'r^un^ooOr— ooroco 

oooooi—LDt— 1— tocvjrovOi— r^ooo^x>csjr^r— mr^oooor^ 
oooooooc— cvjromr^cr»csj'sfr^ocomr^oco^voooo 

OOOOOOOOOOOOOi— •— 1— ouojojcvjrorocororo^ 

oooooododddooooooooddododo 

<4- 

LT) 

II 

2 

<o 

to 

2.147 

2.147 

= 

0.0000 

0.0000 

0.0000 

0.0008 

0.0036 

0.0107 

0.0238 

0.0438 

0.0704 

0.1025 

0.1388 

0.1778 

0.2184 

0.2594 

0.3001 

0.3398 

0.3782 

0.4149 

0.4499 

0.4830 

0.5142 

0.5436 

0.5711 

0.5970 

0.6212 

0.6439 

<4- 

O 

f— 

II 

2 

«ol 

C 

3.220 

3.220 

Ocoo^oo^c^JOOtn^  oro<4*>^^csir^oc\jvoro^cor-o0r—f— 
oocMroir)«?f<oOr-v£>.—  muDCMf^oocou^ 

000»—  ror^CVJCnvOCOr—  OOLOCMOO^OmO^OOCJLOOOr—  ^ 

oooooor— f— c>J^o•e^•^m\i^vor^ooooa^o^o>ooo^— r— 

OOOOOOOOOOOOOOOOOOOOOi— f— F—  r-r- 

4-> 

M- 

IT) 

II 

2 

*o 

2 

o 

c 

CM 

1 

U 

Q) 

i/) 

<o 

6.440 

i 

6.440 

< 

'5* 

■M 

1 

U 

Q) 

(/> 

3: 

oo^co^oocof— voo^r^r^^r- CNi04mroa^»“a>oc\ja>oOr— ^ 
ir)r-^tn^or^0^^oot/)vDcsj^^roro^o^<^a^vooor^roLO 
o«5tcooor^ooooa>vocvj^oa^ooo>r>^^o^Or-voo^oOt— 
oor- c>j^4oa>^c\j^vor-**ooo«“r— 

OOOOOOOf—  r—  r—  I—  CMCMCMCJCsltSJCSJCOOjcvjCNiCMCVJ 

»- 

t: 

o 

•f" 

u 

a 

a. 

(sec) 

ocM^i-'^ooocMTi-^oooocj^voasooj'crvoooocM^ioooo 

r— f— r-p-^CM<MCV4csjcvjrorofnro<n^^^^^mir>mi/>i/ivo 
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L 


1 


. . 2Vfi  sln^kRg 

Bv,  ■ ATn^-iL^ 


"'BV^  ■ 2ALg{l+C^) 

C -_-  ''Vw 

c - J[a_:Li!i^ 

BHj  ALgS^  1 + C^ 


*2  ' f6ALg(l+Cj^) 


-BH3  1 + 


^B^  ~ ^ V kAg  kLg 


= f e-'LB 


4qHk  e~^'^‘-B  sin^kRg 
A kLg  k^A 


B 


sin 


kAg 

4Rb 


All  parameters  in  the  foregoing  equations  have  been  defined  except  C^,  which 
is  the  added  mass  coefficient  in  heave.  The  value  of  is  0.65  and  is  de- 
rived from  the  experimental  data  discussed  in  Chapter  9.  The  coefficients 
which  are  independent  of  the  wave  parameters  are 
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AC 


BV, 


* 0.0758 


ft 


-1 


— ^ = 5.050  X 10-®  L ff  * 

^2  W 


AC 


B5 


= 4.879 


sec 


-2 


AC 


BH, 


= 0.00947  ft'* 


^BH  ” 0.6061 


The  remaining  coefficients  are  given  in  Table  8-4.  The  air  pressure  force  is 
presented  in  Section  8.5  and  the  water  drag  force  in  Section  8.6.  The  analog 
patching  diagram  is  shown  in  Figure  8-3. 


8.5  Air  Pressure  Forces 

The  equations  describing  the  air  pressure  forces  are  developed  in  Chap- 
ter 5 and  the  geometry  of  the  air  passage  is  shown  in  Figure  5-1.  The  re- 
quired scaled  equations  are,  for  - Zg  > 0 

; n\2 


•^aw 

PAL^ 


PAL„ 


I-  • 

z. 


_w 

V 


"B 


aB 


pALj(UC^)  ■ pALj(HC|,) 


w 

V 


and  for  z^  - Zg  < 0 


p 

^aw 

pAi;' 


C V* 
aw 


'w  pAL, 


w 


■^w‘ 

■"b' 

V _ 

V . 
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Table  8-4.  COEFFICIENTS  FOR  BUOY  EQUATION  OF  MOTION 


.0000 
.0000 
.0021 
.0249 
0.0892 
0.1932 
0.3200 
0.4512 
0.5742 
0.6820 
0.7725 
0.8455 
0.9027 
0. 9460 
0.9777 
0.9995 
1.0134 
1.0208 
1.0230 
1.0210 
1.0158 
1.0079 
0.9980 
0.9866 
0.9741 
0.9606 


.6990 
.1311 
-0.6811 
-0.3048 
0.0224 
0.3152 
0.5831 
0.8321 
1.0668 
1.2901 
1.5043 
1.7112 
1.9120 
2.1078 
2.2992 
2.4871 
2.6719 
2.8540 
3.0338 
3.2115 
3.3874 
3.5617 
3.7347 
3.9063 
4.0769 
4.2464 


(sec*') 


0.0662 

0.1682 

0.2780 

0.3745 

0.4507 

0.5065 

0.5449 

0.5696 

0.5836 

0.5898 

0.5902 

0.5864 

0.5797 

0.5709 

0.5606 

0.5495 

0.5378 

0.5258 

0.5138 

0.5017 

0.4899 

0.4782 

0.4668 

0.4558 


H 

(sec  ■') 

“IT 
(sec  ■^) 

0.0466 

0.0000 

0.1738 

0.0000 

0.3676 

0.0003 

0.5782 

0.0083 

0.7864 

0.0455 

0.9219 

0.1278 

1.0368 

0.2490 

1.1174 

0.3892 

1.1699 

0.5282 

1.2005 

0.6519 

1.2146 

0.7534 

1.2163 

0.8309 

1.2091 

0.8858 

1.1955 

0.9210 

1.1773 

0.9396 

1.1561 

0.9452 

1.1329 

0.9405 

1.1084 

0.9280 

1.0834 

0.9098 

1.0580 

0.8875 

1.0328 

0.8624 

1.0079 

0.8355 

0.9834 

0.8076 

0.9596 

0.7793 

0.9364 

0.7509 

0.9138 

0.7229 

r 


PaB  . ^ /] 

z 

w 

pALgd+C^)  pAlgd+C^)  \ 

V _ 

[t\ 

Now,  substituting 


^ = 424.6  in2 
w 

Aj  = 122.7  in^ 
Aj  = 47.17  in^ 
A.  = 2.9  in2 

J 

A„  = 51  in^ 

Ag  = 6786  in^ 
A.  = 25.5  in^ 


Pa-jv.  " 0.002378  slugs/ft^ 

dir 


Wr  " slugs/ft  = 
and,  from  Reference  9,  page  521 


h,,  = 0.34 
h^3  = 0.28 
h3j  = 0.25 
h,3  = 0.56 
h,,  = 0.41 
h^,  = 0.04 

d 1 


the  values  of  the  coefficients  are 


'aw  = 31.89 
e„Caw  = 0.0209 
C^B  =1.980 


EgCaB  = 0.0036 
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The  analog  patching  diagram  Is  shown  In  Figure  8>4. 
8.6  Water  Drag 


The  pressure  drag  Is  discussed  in  Section  4.6.2  and  it  is  noted  that  the 
drag  coefficient  would  be  different  depending  upon  whether  the  buoy  was  rising 
or  falling,  relative  to  the  water.  If  the  relative  velocity  is  defined  as 

V • ii  - ^3 

then  when  this  velocity  is  positive  the  drag  coefficient  would  be  the  forebody 
drag. coefficient  and  when  this  velocity  is  negative  it  would  be  the  base  drag 
coefficient.  Typical  values  are  given  in  Reference  10,  Chapter  3,  and  are 

Forebody  Cp  = 0.90 

Base  Cp  = 0.03 

The  analog  patching  diagrams  to  generate  the  drag  coefficients  and  the  asso- 
ciated drag  forces  are  shown  in  Figure  8-5. 

8.7  Turbine  Torque 

The  torque  produced  by  the  turbine  due  to  the  air  flow  is  discussed  in 
Chapter  6.  Defining 

T^  = torque  scale  factor 

= turbine  angular  velocity  scale  factor 

the  scaled  torque  equation  is 


where  Cj  and  C^  are  given  by  equations  (6-5),  and  (6-7).  Basic  data  for  the 
WATG-4  turbine  is  as  follows: 
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FIGURE  8-4.  AIR  PRESSURE  F 


5 


disk  diameter  20  cm 
blade  span  3 cm 

blade  chord  2.6  cm 

number  of  blades  30 
blade  spacing  2.09  cm 

disk  weight  425  gms 

The  parameters  necessary  for  the  calculation  of  and  are 
Rj  = 3.94  in. 

Ajj  = 1.21  in* 

n = 30 

k = 1 (ref.  8,  pg.  183) 

82  = 25° 
a.  = 15° 

J 

h|j  = 0.95  (assumed) 
and  the  resulting  values  are 

= 0.4285  X 10-'*  lbs-sec*/ft 

Cg  = 0.6220  ft 

The  patching  diagram  for  this  portion  of  the  simulation  Is  shown  In  Figure  8-6. 


8.8  Turbine/Generator  Dynamics 

The  equation  of  motion  for  the  turbine/generator  system  Is  developed  in 
Chapter  7.  Defining  the  additional  scale  factors  as 

= angular  acceleration  scale  factor 
I = generator  current  scale  factor 

the  scaled  equation  of  motion  Is 

[^]  '(^HJ  ■©[%] 
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This  equation  requires  that  the  generator  current  versus  generator  angular 
velocity  be  available  in  the  simulation.  This  relationship  is  simulated  by 
reproducing  Figure  7-1  using  a variable  diode  function  generator.  The  patch- 
ing diagram  for  the  function  generator  is  shown  in  Figure  8-7.  The  generator 
torque  constant,  kg,  is  evaluated  as  follows: 

The  generator  no-load  voltage  characteristic  is 

= kpO) 

0 G 

and  from  Figure  7-1 

Eq  = 69  volts 
u)  = 5000  rpm 


from  which 


kg  = 5^  = .0138  volts/rpm 


or 


kg  = 0.1318  volt-sec/rad 


and  multiplying  and  dividing  by  current 


kg  = 0.1318  watt-sec/ amp 

and  converting  to  foot-pounds 

kg  = 0.1318  (.7376)  = 0.0972  ft-lbs/amp 

The  mass  moment  of  inertia  of  the  rotating  parts  must  be  calculated.  The 
turbine  disk  is  assumed  to  be  a circular  cylinder  of  mass,  m,  and  radius,  Rj, 
so  that 

'^turbine  " 1 

m = 425  grams 
Rj  = 3.94  in 
so 

Jt.rUre  = T (H)  (lb)  (^)’  = 
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The  generator  rotor  is  assumed  to  be  a steel  cylinder  3/4  inches  in  diameter 


and  5 inches  long  so 


and  so 


rotor 


= 1(.283)  J (.75)^(5) 


rotor  " 0-00000948  Ib-ft-sec^ 


~ "^turbine  ^ rotor 
J = 0.001576  Ib-ft-sec^ 


The  generator  damping  constant,  B,  needs  to  be  determined  experimentally  but 
is  estimated  to  be  in  the  range 


0 < B < 0.005  ft-lbs-sec/rad 

The  patching  diagram  for  the  turbine/generator  equation  of  motion  is  shown  in 
Figure  8-8. 


8.9  Parameters  for  a Typical  Simulation 

A typical  WATG  configuration  to  be  analyzed  is  one  with  the  buoy  draft, 
Lg  = 4 feet  and  a water  column  length,  = 15  feet.  These  values  correspond 
to  a buoy  natural  period  of  2.8  seconds  and  a water  column  natural  period  of 
4.3  seconds.  The  coefficients  presented  are  for  the  case  where  the  buoy  is 
forced,  by  a sinusoidal  wave,  at  its  natural  period,  with  a wave  height  of 
H = 1 foot.  The  scale  factors  used  are: 

6 = 10  ft 
V = 20  ft/sec 
A = 50  ft/sec* 

= 500  rad/sec 
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Figure  8-8.  TURBINE/GENERATOR  EQUATION  OF  MOTION 


= 1000  rad/sec^ 

= 20  ft- lbs 
1=5  amps 

The  resulting  potentiometer  settings  are  summarized  in  Table  8-5. 

Various  simulation  runs  were  made  but  these  results  are  not  presented 
since,  until  the  model  can  be  fully  validated,  the  results  would  probably  not 
be  meaningful.  Validation  of  the  model  requires  the  analysis  of  full-scale 
test  data  which  is  not  yet  available.  Chapter  9 discusses  what  data  is  re- 
quired and  presents  the  results  of  an  analysis  of  the  presently  available 
data. 
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9.  VALIDATION  OF  THE  SIMULATION  MODEL 

9.1  Introduction 

In  dealing  with  mechanical  systems,  generally  the  most  troublesome  area 
to  simulate  realistically  is  the  mechanism  of  energy  dissipation.  Usually  the 
dissipation  law  is  not  known  but  assumed,  but  even  if  the  law  is  known,  the 
values  of  the  parameters  involved  are  generally  not  known.  Previous  work  on 
the  WATG  assumed  only  a linear  dissipation  law  but  showed  that  the  power  pro- 
duced is  sensitive  to  the  value  of  the  dissipation  parameter  (damping  ratio  in 
this  case).  The  present  study  assumes  a combination  of  linear  and  second-order 
dissipation  and  estimates  the  values  of  the  parameters  involved  for  both  the 
buoy  and  the  water  column.  The  true  dissipation  mechanisms  and  the  values  of 
the  parameters  involved  needs  to  be  determined. 

The  second  area  involving  validation  is  related  to  the  turbine/generator 
dynamics.  The  torque-speed-air  velocity  relationship  was  derived  in  Chapter  6. 
The  true  functional  relationship  needs  to  be  determined  and  if  the  equation 
derived  is  adequate  the  values  of  the  parameters  Cj  and  Cg  need  to  be  deter- 
mined. Furthermore,  a linear  dissipation  term  and  a linear  electrical  torque 
term  were  assumed  in  the  equation  of  motion  and  this  needs  to  be  validated  and 
quantified. 

9.2  Data  Requirements 
9.2.1  Turbine  tests 

First,  using  the  generator  alone,  a no-load  test  can  be  conducted  wherein 
the  no-load  voltage  is  recorded  as  a function  of  time  as  the  generator  speed 
decays  from  known  initial  values.  Analysis  of  this  data  will  determine  the 
nature  of  the  mechanical  dissipation  torque.  Second,  a similar  series  of  tests. 
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under  load,  together  with  the  known  relationship  between  current  and  RPM,  will 
determine  the  nature  of  the  electrical  dissipation  torque.  Once  these  charac- 
teristics have  been  determined  the  WATG  assembly,  consisting  of  the  generator, 
turbine,  and  air  passage  can  be  mounted  on  a test  stand  such  that  known,  time- 
varying  air  inputs  can  be  introduced.  A means  to  measure  the  torque  delivered 
at  the  generator  rotor  and  to  measure  the  generator  RPM  will  provide  the  data 
necessary  to  deduce  the  torque-speed-air  velocity  relationship. 

9.2.2  Buoy/water  column  tests 

An  instrumented  buoy  is  presently  deployed  and  two  variables  being 
measured  are  buoy  acceleration  and  water  column  displacement  relative  to  the 
buoy.  Analysis  of  the  data  being  gathered  is  complicated  by  the  fact  that  the 
wave  input  is  variable.  If  data  could  be  taken  in  a sheltered  area  of  essen- 
tially calm  water,  a series  of  step  response  tests  could  be  run.  The  buoy 
would  be  lifted  partially  out  of  the  water  and  released.  An  analysis  of  the 
resulting  decay  of  acceleration  and  relative  displacement  would  yield  the 
applicable  dissipation  laws  referred  to  in  Section  9.1.  Some  data  similar  to 
this  has  been  accumulated  and  its  analysis  is  discussed  in  the  next  section. 

9.3  Analysis  of  Available  Data 

9.3.1  Determination  of  the  dissipation  law 

Step  response  tests  using  an  8 x 26  buoy  without  the  WATG  and  with  the 
center  tube  closed  at  the  bottom  were  conducted  by  the  USCG  in  June  and  October 
of  1974.  Acceleration  versus  time  was  recorded  and  five  such  recordings  were 
made  available.  A typical  record  is  shown  in  Figure  9-1.  The  usual  procedure 
to  deduce  a damping  law  from  such  data  is  to  assume  a law  and  then  see  how  well 
(statistically)  the  data  fits  the  law.  If  a linear  law  is  assumed  (dissipative 
force  directly  proportional  to  velocity),  then  a semi-log  plot  of  peak  positive 
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(or  negative)  amplitudes  versus  number  of  cycles  should  yield  a straight  line. 
This  was  the  case  for  all  five  data  records.  The  slope  of  this  semi-log  plot, 
multiplied  by  suitable  constants  yields  the  damping  ratio.  This  value  can  be 
found  by  a least  squares  fit  to  the  line.  The  results  for  the  particular  record 
shown  in  Figure  9-1  yielded  a value  of  0.0461  when  the  positive  peaks  were  used 
and  a value  of  0.0463  when  the  negative  peaks  were  used.  For  a linear  system 
the  same  value  of  damping  ratio  should  be  obtained  if  it  is  computed  by  taking 
the  ratio  of  a positive  peak  and  the  next  negative  peak  (disregarding  the  sign) 
or  a negative  peak  to  the  next  positive  peak.  When  this  was  done,  it  was  found 
that  the  values  obtained  were  not  the  same.  The  implication  of  this  is  that 
the  energy  dissipated  is  different  depending  upon  whether  the  buoy  was  moving 
out  of  the  water  or  into  the  water.  A complete  analysis  using  all  five  sets 
of  data  yielded  the  following  values; 

Moving  out  of  the  water,  damping  ratio  = 0.0664 

Moving  into  the  water,  damping  ratio  = 0.0319 

Using  these  values  in  the  linear  model,  peak  amplitudes  were  computed 
and  compared  to  the  measured  values.  A regression  analysis  yielded  a coeffi- 
cient of  determination  greater  than  0.99  which  means  that  less  than  1%  of  the 
difference  between  the  calculated  and  measured  data  is  not  accounted  for  by 
the  assumed  law. 

The  physical  explanation  for  this  behavior  would  appear  to  be  that  when 
the  buoy  is  moving  into  the  water,  viscous  forces  are  generated  but  when  moving 
out  of  the  water,  a wake  is  developed  which  generates  additional  viscous  forces. 
This  form  of  energy  dissipation  might  be  termed  bilinear  viscous  damping;  the 
dissipative  force  being  directly  proportional  to  velocity  but  with  the  constant 
of  proportionality  having  one  value  for  positive  velocity  and  another  value  for 
negative  velocity.  Two  questions  remain  with  reqerd  to  the  present  simulation; 
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first,  is  this  law  applicable  to  the  buoy  when  the  center  pipe  is  not  capped 
and  second,  what  is  the  dissipation  law  for  the  relative  water  motion  in  the 
center  pipe.  These  questions  can  be  answered  by  conducting  the  tests  outlined 
in  Section  9.2.2. 

9.3.2  Determination  of  the  added  mass  coefficient 

The  step  response  data  yields  a damped  natural  period  of  the  buoy  tested 
of  2.7  seconds.  Since,  as  was  shown,  the  damping  is  small,  this  can  be  taken 
as  the  undamped  natural  period.  The  buoy  natural  frequency  is,  therefore, 

ug  = ^ = = 2.33  RAD/SEC 

Now,  since 


and,  for  the  buoy  tested  Lg  = 3.61  Feet 
* 0.65 

which  is  the  value  used  to  calculate  the  coefficients  presented  in  Chapter  8. 
Since,  for  the  buoy  tested,  the  center  pipe  was  capped,  future  tests  as  pro- 
posed in  Section  9.2.2  will  determine  any  necessary  correction  to  this  value. 

9.3.3  Modification  of  the  Simulation  equations 

The  validation  procedure  may  require  that  certain  forces  in  the  basic 
equations  be  modified  or  that  the  numerical  value  of  coefficients  used  be 
modified.  The  forces  involved  are: 

*^Wg  = water  column  friction  force.  Section  3.3 
» buoy  friction  force.  Section  4.6.1 
Dp  = buoy  pressure  force.  Section  4.6.2 
Dg  ■ buoy  radial  dissipation  force.  Section  4.6.3 


69 


T 


3 


turbine  torque.  Section  6.2 


Bo)  = turbine/generator  mechanical  dissipation  torque.  Section  7,2 
kgi  = turbine/ genera tor  electrical  dissipation  torque.  Section  7.2 

If  the  damping  law  discussed  in  Section  9.3.1  proves  applicable,  the 
simulation  equations  would  be  significantly  less  complex. 
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10.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  objectives  of  this  project,  as  stated  in  the  Introduction,  have  been 
met.  The  buoy,  water  column,  air  passage,  turbine  and  generator  have  been 
modeled  and  analog  computer  patching  diagrams  presented.  The  available  full 
scale  data  has  been  analyzed  and  this  analysis  indicates  that  modifications 
in  the  model  will  probably  be  required.  This  conclusion  cannot  be  made  more 
definite  since  the  data  analyzed  was  for  the  buoy  alone  without  a water  column. 
It  is  likely,  however,  that  even  with  a water  column,  the  dissipation  law  will 
be  proved  to  be  bilinear  thus  simplifying  the  equations  as  presented  herein. 

It  is  recomnended  that  any  continuing  work  on  this  project  include  the 
following: 

(1)  The  tests  of  the  turbine  and  generator  outlined  in  Section  9.2.1. 

(2)  The  buoy /water  column  tests  as  outlined  in  Section  9.2.2. 

(3)  Modification  of  the  simulation  equations  as  required  by  the  afore- 
mentioned tests. 

(4)  Generation  of  Response  Amplitude  Operators  from  the  simulation  model. 

(5)  Comparison  of  the  RAO's  generated  by  the  simulation  with  those  to 
be  available  from  the  deployed  buoys  now  under  test. 

(6)  Presentation  of  the  results  in  a form  suitable  for  the  USCG  digital 
Design  Integration  Model  currently  under  development. 

The  possibility  exists  that  analysis  of  the  pure  heaving  motion  may  not 
be  adequate  to  describe  the  actual  buoy  behavior.  If  this  should  be  the  case, 
then  the  model  would  have  to  be  modified  to  include  the  coupled  tilt-heave 
dynamics. 
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